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Abstract 



We propose using decimation to estimate vehicle speed 
for adaptive channel estimation in pilot-aided transmission 
systems like DS-CDMA ones. It is shown that by passing 
the received pilot signal through bandpass filter after deci- 
mation and computing its output energy, current velocity of 
mobile station can be accurately estimated. In addition, the 
proposed speed estimation unit is simple to implement. This 
allows us to use it accompanying computationally efficient 
pilot filter for forward link receiver, a structure for adaptive 
channel estimation based on the proposed velocity estima- 
tion unit is addressed. The optimal coefficients of pilot filter 
is selected from a precalculated bank of coefficients based 
on the estimated mobile speed zone. Procedures for deriv- 
ing the optimal coefficients and thresholds for the speed es- 
timation unit are developed. The proposed method provides 
up to 0. 75dB improvement compared to conventional non- 
adaptive ones. 



1. Introduction 

In the mobile communication systems employing pilot 
tone as a phase and amplitude reference for demodulation 
on fading channels, the pilot filter is used to extract the ref- 
erence. For example, a bank of lowpass filters is employed 
in DS-CDMA system to estimate the channel response for 
each finger prior to rake combining. The receiver extracts 
the pilot using pilot filters with fixed bandwidth. In the for- 
ward link, some computationally efficient pilot filters like 
a moving average filter and a first-order EIR filter are pre- 
ferred for simple implementation. The bandwidth of pi- 
lot filter should be wide enough to accommodate the fade 
spectrum (the maximum Doppler frequency fo). If it is too 
narrow, the estimated channel response cannot trace the ac- 
tual channel fluctuation which leads to the appearance of 
an error floor. In contrast, if it is too wide excessive noise 
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appears in the output of pilot filter. The optimal pilot filter 
for a given fry can be obtained using Wiener filter theory. 
However, the Wiener filter for a fixed fry cannot be optimal 
in practice because fry varies as the mobile speed changes. 

The performance penalty caused by fixed bandwidth can 
be reduced by introducing adaptivity to pilot filtering so that 
the bandwidth of pilot filter is adjusted to changes in the 
Doppler. Some adaptive filtering schemes were previous- 
ly proposed [l]-[2]. Use of the LMS algorithms based on 
multi-dimensional adaptation [1] is too complex to be used 
in the forward link. More reliable adaptation algorithm pro- 
posed in [2]. When the magnitude response of a pilot filter 
is not flat in passband like a moving average filter and a first- 
order IER filter, it is not adequate. Many literatures have 
shown to estimate velocity of mobile and corresponding fry 
[3]- (7]. These methods are computationally inefficient or 
unreliable to be used in practical systems. For example, the 
performance of the method based on level crossing detector 
[7] is severely degraded under AWGN. 

Another performance penalty of the currently used chan- 
nel estimation method in DS-CDMA systems is caused by 
using the prediction-type channel estimator — group delays 
caused by pilot filtering is not compensated in traffic chan- 
nel decoding. Compensation of this kind of group delays 
degrades the performance of power control due to the addi- 
tional delays in its loop. 

In this paper, a new mobile speed estimation method us- 
ing decimation is introduced. This estimation unit is re- 
liable and simple to implement. Using the proposed mo- 
bile speed estimauon method, an adaptive channel estima- 
tion scheme reducing or eliminating the performance penal- 
ties mentioned above is addressed. Up to 0.75 dB gain was 
gained using the proposed scheme. 

2. System model 

In this paper, we only consider the pilot channel for chan- 
nel estimation in DS-CDMA system. The performance of 
channel estimauon could be improved using not only pilot 
channel signal but also traffic channel signal. Then trans- 



0-7803-4507-0/00/$ 10.00 ©2000 IEEE 



2839 



VTC 2000 



5DOCID: <XP 1033051 A_J_> 



mitted signal s(t) assumed to be QPSK modulated and QP- 
SK spreaded is given by 



5(0 



(1) 



+ 3 



{ai(t) + a Q (t)}] 



+ traffic channel signals 

where, E CyP is pilot signal energy per pseudo noise (PN) 
chip sequence, and a/(£) and a Q (t) is I-channel and Q- 
channel PN sequence, respectively. In Eq. (1), transmitter 
filter h(t) is absorbed in aj(t) and a Q {t). To simplify analy- 
sis, we shall initially assume frequency non-selective fading 
and single user environment, then received pilot signal r(t) 
can be written as 



r(t) = \R h (T)a^^{ar(t)-a Q {t)}cos<l> (2) 



+ 



\Rh{r)a^^{aj{t) + a Q (t)} sin 4> 

+ l^h(r)a^^{ ai (t) + a Q (t)} cos<£] 
+ terms bearing traffic channel signals 
+ \ni{t) +jl nQ {t) 

where, a and <p is amplitude and phase response of Rayleigh 
fading, respectively, R h (r) is autocorrelation function of 
h(t) under the assumption of using matched filter at receiv- 
er, and r is timing error. n T (t) and UQ{i) are white Gaus- 
sian noise, and their single sided spectral densities are A' 0 's. 
Under the assumption that there is no timing error (r = 0), 
the terms bearing traffic channel signals in Eq. (2) can be 
discarded. Moreover, ^(0) = 1. Finally, pilot despreaded 
received signal y p (n) prior to pilot filtering is expressed as 



y p (n) = a(n)N^/2E CiP cos 4>{n) 
+ joL{n)Ny/2E CiP sin <f>(n) 
+ n i( n ) + jn' Q (n) 



(3) 



where, N is a spreading factor and E[n'j(n) 2 ] = 
E[riQ(n) 2 ] = NN Q . The channel estimation is performed 
by manipulating y p (n) given in Eq. (3) on symbol rate. 

3. Proposed architecture 

The proposed adaptive pilot filtering scheme is based on 
the estimation of mobile speed. Tap coefficients of pilot fil- 
ter are updated to adjust its bandwidth to accommodate f D 



corresponding to the estimated speed. Determination of the 
type and order of pilot filter p(z) depends on applications. 
However, it should be simple to implement for forward link 
receiver. Throughout this paper, we assume that the follow- 
ing first-order IIR filter is used for pilot filtering — it is more 
relevant to use a first-order IIR filter for adaptation than to 
use a moving average filter 



p{z) = 



1 - (1 -a)z~ l 



(4) 



where, a is a power-of-two coefficient so that it is simple to 
implement. Then, it is the only a to be adaptiveiy updated. 
In addition, the speed of mobile stations doesn't need to be 
estimated exactly because only a few discrete values of a 
is available. The mobile speed to be estimated is divided 
into a few number of target speed zones, and a set of the 
optimal coefficient a for each speed zone is precalculated 
and stored. In results, each speed zone has a different pilot 
filter with a different bandwidth (coefficient a) optimized 
for each velocity. The best one is automatically selected 
according to the estimated speed zone. The proposed archi- 
tecture is shown in Fig. 1 . The group delay dt caused by 
pilot filtering is compensated for traffic channel decoding. 
In other hands, a separate path without this compensation is 
provided for power control. In general, it is impossible to 
determine the exact value of dt because IIR filter has non- 
constant group delay. In this paper, the group delay at DC 
component is selected for dt . The performance degradation 
caused by this approximation is negligible. It can be easily 
shown that^ = [1/aJ — 1. 

3.1. Estimation of mobile speed 

The proposed mobile speed estimation method shown in 
Fig. 2 is based on searching fr> by using decimation pro- 
cessing. As mentioned above, only a few number of speed 
zones are sufficient to accommodate a finite number of dis- 
crete value of a. Each branch with a M-to- 1 decimator in 
Fig. 2 is corresponding to one mobile speed zone. The alias- 
ing caused by decimation is reduced by passing output of 
pilot filter to M-to-1 decimator. By decimating the pilot 
signal by M, its effective bandwidth is widened by a fac- 
tor of M. In results, there exists some M"s which make the 
averaged power of bandpass filtered pilot signal eQ(n) larg- 
er than the predetermined threshold U as shown in Fig. 2, 
where integer i represents each branch. The bandpass filter 
with high Q-factor here is a simple device for energy mon- 
itoring. The current speed zone of mobile station is the one 
corresponding to the minimum value of M Accordingly, 
the decimation in Fig. 2 should be continued until eQ(n) 
becomes larger than U at i th branch for the first time. If 
there is no M making ep(n) larger than any t^'s, the cur- 
rent speed zone is determined to the slowest velocity. In 
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facu one more threshold t' { exists at each branch to ensure 
backward estimation. That is, the current speed zone moves 
up to faster speed zone of (i — l) th branch if eQ(n) > t\. 
The details for deriving the optimal threshold will be ad- 
dressed later. 

To illustrate the proposed mobile speed estimation u- 
nit, consider the following DS-CDMA example: chip rate 
is 3x1.2288MHz, data rate is 9600bps, symbol rate is 
14400sps, processing gain is 256, and carrier frequency is 
assumed at 2GHz. The center frequency of the bandpass 
filter in Fig. 2 is fixed at 0.1. Let the current speed of mo- 
bile be 50 km/h y then the maximum Doppler frequency 
fo is 92.59Hz which is corresponding to fo = 0.0064 
in normalized frequency. After 8-to-l decimation, fo be- 
comes 0.0514. Because fo « 0.1, eg(n) << to- There- 
fore, decimation should be continued. After the next 10- 
to-I decimation (M = 10), f D = 0.0643 which leads 
to ec?(n) << t y again. After decimating by 16 totally, 
f D = 0.1029 and cq{ti) becomes larger than t 2 for the 
first time. So, the estimated current velocity is between 30 
km/h and 60 km/h. As velocity increases to 70 km/h, 
eq(n) becomes larger than t 2 , so that current velocity is set 
to the desired speed zone (i = 1). In the same way, the 
last zone is selected as the speed decreases to 20 km/h as 
e Q( n ) < *2- Choosing the number of decimation branches 
and a decimation ratio in each branch depends on applica- 
tions. 

As total power of received signal increases, eQ(n) can 
be also indefinitely increased, so that thresholds seem to be 
meaningless any more. In practice, however, it does not 
happen because automatic gain control (AGC) is employed 
in every cellular system where dynamic range of received 
signal is about 80dB. Otherwise, it can be also avoided by 
normalization. 

3.2. The optimal pilot filter 



y p (n) in Eq. (3) can be rewritten as 
1 



2/» = 



1 



(5) 



=n(n) 



where, /3(n) = a(n) cos <j>{n) + jcr(n) sin<£(n), cr^ is the 
variance of f3(n), and n(n) = n' 7 (n) + j^Q(n). The spec- 
trum of y$(n) is given by 



S 0 (f) = * f2 
and its autocorrelation function is 



(6) 



(7) 



Jo( ) is Bessel function of the first kind with zero order. The 
optima] coefficient a of pilot filter p(n) is obtained by min- 
imizing the following one dimensional optimization prob- 
lem with respect to a. 



F(a) 

= E\\p(n) * *y' p (n) - 0(n - d t )\\ 2 
* ^p(k)p(l)Jo{27rf D (k - l)T) 

k I 

N 0 



(8) 



+Jo(0) + 



NE C 



-2^p(i)A(2T/D(t-4)T) 
k 

where, T is the symbol duration and d t is the group de- 
lay explained in Section 3. This optimization problem can 
be easily solved numerically because only a few values of 
a are available (a is a power-of-two coefficient). For in- 
stance, if 2 _1 , 2~ 2 , and 2~~ 3 are available for the value of 
a, the optimal one for given fo can be obtained by evalu- 
ating and comparing the value of F(a) each other. Notice 
that the optimal coefficient a is a function of pilot signal-to- 
interference ratio (PSIR) E CyP a^/N Q as it is expected. 

3.3. Determination of threshold ti and t\ 



The pilot filtered y p (n) is 



y p = Ny/2E CtP p(n) * *p(n) + n(n) 



(9) 



where, E[n 2 (n)] = 2NN Q £ fc pW 2 - Consider i th branch 
with M-to-1 decimator. This branch is corresponding to the 
mobile speed of between v\ and v!>. After M-to-1 decima- 
tion and bandpass filtering, y p becomes 



MN^y2E^,0{Mn) * p(Mn) * b(n) + n(n) (10) 

where, b(n) is the impulse response of bandpass filter in 
Fig. 2. We assume that the magnitude response of p(M n) 
in Eq. (10) is almost flat in passband [0, fo] and n is white 
for simplicity. Then, the total output energy is given by 

*qUd) » 2M 2 N 2 E CtP ay (11) 
E Hk)b(l)M2nMf D (k - 1)) 

A: / 

+2MNN 0 Y,P(V 2 Y, b W 2 
k k 

In results, U = eQ(f l D ) and % = ep(/£), respectively, 
and f l D = v\/X (/£> = v+/X). AGC or normalization can 
be used for preventing eQ(fo) to increase so large as men- 
tioned Section 3.1. 
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3.4. Practical considerations 

Despite all a, U's y and t\ are functions of PSIR, they are 
quite robust to it. It seems to need a huge table of them for 
every possible PSIR, but only a few number of PSLR's are 
sufficient because a few number of possible a's are avail- 
able. Roughly estimated PSIR is also accessible in the re- 
ceiver. 

4. Simulation results 

To demonstrate the efficiency of the proposed method, 
the proposed method was applied to the forward link of the 
DS-CDMA system considered in Section 3.2. Convolution- 
al code with the rate of 1/3 and the constraint length of 9 
is employed. Fast power control was applied with the rate 
of 800 Hz and power control error of 10%. 20% of the 
cell power is assigned to the pilot. Jakes* model was used 
for generating frequency selective Rayleigh fading channel. 
The first order IIR pilot filter is employed for pilot filtering 
and the coefficient a is updated frame by frame. The mobile 
speed is divided into five speed zones of {[0, 30 km/h), 
[30, 60 km/h), [60, 90 km/h), [90, 120 km/h), and 
[> 120 km/h)} as shown in Fig. 2. Corresponding stored 
coefficients a are {16,8,4,4,2} with E CyP /N Q = -13dB. 
Notice that a's in 3 rd and 4 th zone are same. We could 
merge two zones but we don't for the proper operation in 
the different PSIR's. 

In general, frame error rate less than 2% is acceptable for 
voice decoding. Gains in the averaged transmitted power to 
achieve the frame error rate less than 2% over conventional 
one with a fixed a = 1/4, are shown in Fig. 3 for several 
vehicle speed in km/h. Up to 0.75dB gain (yielding up to 
18.7% of capacity increase) is obtained by introducing the 
adaptivity in channel estimation. The gain drastically in- 
creases as the number of multi-path or the allowable com- 
plexity in pilot filter increases. The worst case convergence 
time is occurred when the velocity is between 0 km/h and 
30 km/h, because it needs more than two frames to get 
a reliable eQ(n) — as decimation ratio M increases, a few 
number of symbols in a frame are available for computing 
it. However, the performance degradation is negligible be- 
cause the amount of variation in vehicle speed for two or 
three frames is very small even if an acceleration of 0-100 
km/h in 5 seconds. The details of analysis is not addressed 
here because of space limitation. 

5. Conclusion 

An adaptive channel estimation method relying on esti- 
mating the velocity of mobile station for pilot-aided trans- 
mission systems was developed. It was shown that the sig- 



nificant improvement in system performance could be pro- 
vided by introducing adaptivity in pilot filtering. Moreover, 
the proposed architecture was simple for implementation, 
so that it was appropriate for forward-link receiver. 
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Figure 1. Proposed architecture for channel 
estimation in DS-CDMA systems. 
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Figure 2. Proposed mobile speed estimation 
unit. 
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